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, We find that tunnel junctions are not needed to provide single electron effects in a metallic island. 

' Eventually the tunnel junction may be replaced by an arbitrary scatterer. It is important that even 

a diffusive scatterer provides a sufficient isolation for single electron effects to persist. To formulate 
this in exact terms, we derive and analyze the effective action that describes an arbitrary scatterer. 
We also consider the fluctuations of the effective charging energy. 
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It is well-known that electric charge of an isolated piece of conducting material can only take discrete values 
corresponding to the integer number of electrons in there. This property persists if this isolated piece, the island, is 
connected to electron reservoir by means of a resistive tunnel junction It is the recognition of this mere fact that lead 
r ^ . to outburst of the entire field of single electron phenomena. 

The single electron effects are best visible provided the conductivity of the tunnel junction is much smaller than 
^ . the conductance quantum Gq = e 2 /2nTi. The ground state energy as a function of induced charge q is given by 
minimization of Coulomb energy, Ec(n + q) 2 , with respect to discrete charge n. The result is periodic in q with a 
period e. Q The analysis of the reverse case, G ^> Gq, requires advanced theoretical methods. Despite the partial 
controversy in results, all authors agree that in this case the ground state energy retain the periodic q dependence, 
that manifests the Coulomb blockade. The effective charging energy, Ec, that is, the g-dependent part of the ground 
state energy, is suppressed by a factor of exp(— G/2Gq) in comparison with Ec- 

Still the analysis has been restricted to tunnel junctions. Next step has been made in |^|| where Coulomb blockade 
has been studied in the situation where the isolation is provided by a quantum point contact with almost perfect 
q i transparency. It has been shown that the charge quantization survives. Albeit the charging energy is strongly 
Q suppressed vanishing to zero at perfect transmission. 

In this paper, we construct a general theory of Coulomb blockade that can embrace tunnel junctions, quantum 
point contacts, diffusive conductors and eventually any type of scattering. 

The results are as follows. Charging energy vanishes only for perfect point contacts. For very wide class of conductors 
that have conductivity G ^> Gq, the charging energy is exponentially suppressed, ln(E c / E c ) oc —o:G/Gq, a being 
dimensionless coefficient depending on the type of the conductor. For disordered conductors, for instance, diffusive 
ones, the charging energy strongly fluctuates. This happens even if the fluctuations of the conductance are small. 

The most equivalent mathematical framework to describe the charging effects in full has been proposed by Schon 
and Zaikin. JtJ They present the partition function of the system in the form of the path integral over the field 
OV <p{r){[3 = n/T), 
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The form of the last term presumes that the partition function can be presented as a sum over topological sectors 
O ■ that are labelled by an integer W — dTip(r) winding number of fir), 

Z = J2 Z wexp(-i2TrWq/e) (2) 

. Here Zw does not depend on q. 

The action C sc describes tunnel junction. It can be evaluated with using Tunneling Hamiltonian method that gives 
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This form shows that the tunnel junction is quite different form a linear resistor which is described by a form 

r — — Fa Fa ' ^ Z ^ T '^ m 
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which is bilinear in ip. It has been frequently assumed that a coherent diffusive conductor can be described by (Q) 
and consequently exhibits no charging effects. We show below that it is not so. However, the relation (^) holds for 
an arbitrary conductor in the limit of small ip. 

We sketch the derivation of C sc for an arbitrary conductor. Our basic assumptions are: (i) no inelastic scattering 
occurs in the conductor, (ii) the conductor is sufficiently short, EcT trav /h -C 1, T tra v being typical traversal time 
throught the conductor, (iii) the island is sufficiently large, so that Eq greatly exceeds the average level spacing in 
the island. Two first assumptions allow us to characterize the conductor by an elastic scattering matrix disregarding 
retardation effects. The third assumption allows us to disregard coherence between the electrons transmitted to and 
coming from the island, so that it can be regarded as an electron reservoir. 

These assumptions correspond to a fermionic action of the following form: 

-£= dr(^2 / dxxi{x,r) [d T + iv n d x ] Xn(x , r) + i/jI(x,t) [d T - iv n d x ]ip n (x, r) (5) 

J0 n =i-J-°° 

+ d,xxi{x,r) [d T - iv n d x ]xn(x,T) +^(x,r) [d T + iv n d x ]ip n (x,T) 

+E c (Q(r)-q/e) 2 ) 

Here the island is on the right (x > 0), n and m label transport channels in the island and the reservoir respectively. 
X stands for fermion fields coming to the scatterer and ip stands for outgoing modes. The scatterer is completely 
characterized by the scattering matrix Ski that set a boundary condition for ip and x 

^fc(0) = £3wxi(0). (6) 
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Here k, I label modes on both sides of the scatterer. 
The charge in the island is given by 

Q(t) = ^2 / dx (xL0> T )Xm(2:,"r) +ip} n (x,T)ip m {x,T)) (7) 

m=l J ° 

To proceed, we perform a Hubbard-Stratanovitch transform on interaction term introducing a new variable (p(r), 
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Ec(Q(r) - q/ef) - -^L + zQ(r)0(r) - iq/e^r) (8) 

The resulting action is quadratic in fermions so that they can be integrated out and the action can be represented 
as a functional of (p(r). It is constructive to do this calculating Green functions of the fermions in the presence of 
field (p and scattering potential. Special attention shall be given to the fact that in one dimension the Green function 
G(x, r, x',t'; [<p(t)]) is not continuous at coinciding arguments. So that the problem shall be regularized by letting 
scattering to occur in small but finite region of space. We present the action as a series in AG = Gq(x, r, x', r'; [p(t) = 
0]) — Gq(x, t, x' , t'; [<£>(t)]) where Go refers to Green functions of disconnected island and reservoir, since AG is 
continuous. Then the action can be reduced to the trace of logarithm of an operator, 

- C sc = Tr n , r m(l - (I - f)§ + cxp(-iip)§~ exp(i^) — (I — r + )6~ exp(-iip)§ + exp(iip)) (9) 

The operator is defined in a Gilbert space which is the direct product of the Gilbert space of the transport channels 
in the island and the Gilbert space of fermion functions in imaginary time. The operator f is the reflection matrix 
in the island. The operators fl* are projection operators on positive (negative) Matsubara frequencies, 6^ = 9(±e n ). 
Although 8 ± obviously commute with f, f + , they do not commute with multiplication by exp(±z</?(r)). This generates 
a complicated operator algebra that makes a complicated problem even the evaluation of C sc . 

The expression (^|) can be explicitly evaluated in two limits: f — > I and <p — > 1. There, we successfully reproduce 
Eq. |^ for tunnel junctions and Eq. ^ for an arbitrary scatterer in linear regime. 

Fortunately, a very important part of analysis of the action (^) can be done exactly. We are able to find the 
minima of (|p) in each topological sector and thus give a quantitative estimate of effective Coulomb energy in the limit 
G 3> Gq. We consider the configurations of ip(r) of the following form: 
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ex P M = [] T — -L. (10) 
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Here u = exp(i27rr//3), Zi are complex parameters. Zj can be viewed as coordinates of AT (anti)solitons in the plane 
of complex u (see Fig. 2). If |zj| — > 1, these configurations correspond to sets of Korshunov's solitons. |l We are 
interested in configurations where all solitons are of the same sign. In this case either |zj| < 1 for all i or \z%\ > 1 so 
that exp(iip) is an analytical function of u either within or beyond the unitary circle. The winding number W — ±N. 
Using the methods of analytical function theory we show that these configurations indeed minimize the action in the 
corresponding topological sector. The minimum does not depend on Zj and equals 



- In det (ff+) | W | + - 



. w = - In det (ff+) | W | + ^ lndet(f/f+)W. (11) 



The second term is imaginary and can be viewed as a trivial shift of induced charge q: q — ► q + ilndet(f/r+). The 
first term is of importance since it describes the suppression of statistical weight of topogical sectors with W ^ in 
comparison with the trivial sector. It has been shown in |||| that the suppression of these statistical weights leads to 
suppression of effective charging energy. This allows us to write down a simple formula for effective charging energy 

Ec^EcWrH 2 (12) 

n 

where R n are eigenvalues of the reflection matrix ff + . This formula is valid provided the suppression is big. Similar 
relation has been obtained by Flensberg (g) in a much more restrictive framework. In the limit of almost perfect 
transmission, R — > 0, we reproduce the results of Matveev ||. 

Recent theoretical advantages allow to characterize R n of a scatterer/conductor of virtually any type (see || for 
review). This makes the relation ( [l2|) easy to use for concrete examples. From now on, we will concentrate on diffusive 
conductor in the limit G Gq. It is a disordered conductor, so that it is characterized by distribution of i?„, or 
transmissions T„ = 1 — R n . It has been shown in fllOfl that the transmission distribution of a diffusive conductor 
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depends only on its conductance, p(T) = Gj (2GqT\J1 — T). We average logarithm of ( |12|) with this distribution to 
obtain 

Ec/EcxeM-^)- (13) 

This is the main result of this work. The diffusive scatterer of the same resistance as a tunnel junction suppresses 
Coulomb blockade much more efficiently. To give some numbers, let us choose 1 /G = 4 kOhm. In this case, suppression 
factor is about 25 for a tunnel junction and almost 3000 for a diffusive conductor. 

Below we consider the fluctuations of Eq and the effect of weak localization. To make a qualitative estimation, we 
note that the fluctuation of G is of the order of Gq. The weak localization correction is of the same scale. Therefore, 
the fluctuations of an exponential like ( |l3| ) must be of the order of its average value. The same should hold for the 
weak localization effect. It is remarkable that quantitative consideration gives even bigger values. 

This quantitative treatment can be performed along the lines of |llj and |L2]]. There are formulas that can be 
directly applied to the quantity of interest \n(Ec / Ec) — Sn m ^™- ^ appears that both the fluctuation and the 
localization effect are dominated by a contribution of the universal cooperon-diffusion mode, the one which provides 
Wigner-Dyson statistics of closely spaced transmission eigenvalues, jllj The contribution of this mode logarithmically 
diverges at very small R and shall be cut off at R ~ Gq/G, average value of transmission spacing. For pure statistical 
ensembles, the fluctuation is given by 

« \n 2 {E c /E c ) »= HG/G Q ) (14) 

where Nd c is the number of massless cooperon and diffusion modes. It ranges from 1 to 8. The weak localization 
correction is 

< HE C /Ec) > wi = -^f HG/G Q ) (15) 

where N w i is 2, 0, -1 for simplectic, unitary and orthogonal ensemble respectively. 

Experimentally, the fluctuation and weak localization effect are identified by their magnetic field dependence. 
Following [pT| we introduce dimensionless parameters i]H,r]so to characterize magnetic field and spin-orbit interaction. 
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We disregard influence of magnetic field on spin. Correlator of two Eq taken at different values of magnetic field 
reads 

« HE c { m )/E c ) HE C ( V ' H )/E C ) »= -\h\tjI - V 'h 2 \ 4 ((vh + r,' H f + 4 f{{ m - v'h? + 4of)- (16) 

This equation is valid provided G/Gq <C \t]h — 1 i'h\>''1h,v'h>' i 1so *C 1. It shows that the actual value of E c can be 
changed by an order of magnitude by fairly small change of magnetic field 5t]h — Gq/G. At a bigger Srj H effective 
charging energy exhibits power-law correlations. 

Magnetic field dependence of the weak localization correction is given by 

< HE C /Ec) > wl = \$HnH + vlo) -Hv%))- (17) 

To conclude, we have shown that the isolation required for discrete charge effects can be provided by any constriction 
which is not ideally ballistic. We have discussed suppression of the effective charging energy by a diffusive scatterer 
and found gigantic fluctuations of this quantity. 

I am indebted to G. W. E. Bauer, P. W. Brouwer, S. E. Korshunov, Y. Gefen and many others for very instructive 
discussions of the results. This work is a part of the research programme of the "Stichting voor Fundamcnteel 
Onderzoek der Materie" (FOM) and I acknowledge the financial support from the "Nederlandse Organisatie voor 
Wetenschappelijk Onderzoek" (NWO). 
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